Modulated phases in magnetic models frustrated by long-range interactions 
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We study an Ising model in one dimension with short range ferromagnetic and long range (power 
law) antiferromagnetic interactions. We show that the zero temperature phase diagram in a (lon- 
gitudinal) field H involves a sequence of up and down domains whose size varies continuously with 
H, between — H c and H c which represent the edge of the ferromagnetic up and down phases. The 
implications of long range interaction in many body systems are discussed. 



I. INTRODUCTION 

It is known that long range interactions lead to signif- 
icant changes in the behavior of interacting many body 
systems. Thus, though the Ising model in one dimension 
with a short range ferromagnetic interaction does not ex- 
hibit a phase transition, the corresponding model with a 
ferromagnetic interaction that falls off as a power law 
1 /r x has a phase transition at non-zero temperature^ for 
A < 2. A similar effect is found for the Ising spin glass in 
one dimension, which is unfrustrated with nearest neigh- 
bor interactions, but becomes sufficiently frustrated with 
a power law interaction so that for A < 1 a finite T tran- 
sition is attained 2 -^. 

For short range models with a finite transition tem- 
perature, addition of long range interactions having a 
power-law fall off can lead to changes in the universal- 
ity class of the phase transition. For sufficiently small 
power law exponent, critical exponents are found to vary 
continuously with the power law exponent^. 

There have been several studies of the Ising model in 
higher dimensions, frustrated by Coulomb^ or dipolar^ 
long range interactions, but without a magnetic field. 
Here we examine the effects of a frustrating long range 
interaction on the phase diagram of a ID Ising model 
at zero temperature in the presence of a magnetic field. 
This work is motivated by the proposal by Spivak and 
Kivelson^ (and generalized with JameiiS) that the pu- 
tative first order phase transition between the Wigner 
crystal and Fermi liquid phases of the interacting elec- 
tron gas in two dimensions at T = is pre-empted, 
due to the long-range Coulomb interaction, by a series 
of "microemulsion" phases characterized by phase sepa- 
ration on a mesoscopic scale. In general, a system with 
a long range interaction that frustrates the order favored 
by a short range interaction will not macroscopically sep- 
arate into the phases of the unfrustrated system, once the 
long range interaction is strong enough. This is because 
the short range interaction energy increase due to having 
mesoscopic domains is smaller than the long range in- 
teraction energy increase due to having macroscopic do- 
mains. (Thus, in such systems, the Maxwell construction 
for determining phase separation must be generalized^.) 

In this work, we demonstrate that a Coulomb frus- 



trated Ising model in a spatial continuum in one dimen- 
sion (and its generalization to other power laws) admits 
analytical solution at zero temperature. We find that this 
system possesses a regime exhibiting modulated phases 
[i.e. mesoscopic phase separation), with a period that 
varies continuously with applied magnetic field. 



II. THE MODEL 

We study a one-dimensional frustrated Ising model 
given by adding to the familiar ferromagnetic Ising chain 
in a magnetic field, where only nearest neighbors inter- 
act, a competing long-range antiferromagnetic interac- 
tion with a power law fall-off. We investigate the model 
in a one-dimensional continuum, where its Hamiltonian 
is given by: 



U H = J / dr 



da(r) 



dr 



oo 

H J dra(r) 



Q 



+ — / dndrj v(r)a(ri)a(rj) . (I) 



Here o~(r) = ±1 is a bi- valued function of the posi- 
tion r (Ising spin), J and Q are parameters representing 
the strengths of the short-range ferromagnetic and long- 
range antiferromagnetic interactions respectively, and H 
is the strength of a uniform magnetic field. In this paper, 
we take the long-range antiferromagnetic interaction as 
v(r) = l/(r + a) A , with exponent A > 0, and an ultravio- 
let cutoff a that must sometimes be retained to eliminate 
divergences. In the case H = 0, this model has been 
solved by Giuliani et a/.—, and for H = and v(r) equal 
to the inverse Fourier transform of the inverse Laplacian, 
Grousson et al*£ have analytically solved the model in 
three dimensions. 

It proves helpful to perform a Legendre transformation 
on the energy and work at fixed average spin density a 
instead of fixed magnetic field strength H, especially for 
A < I where the interaction energy density is infrared di- 
vergent for a ^ 0. Since the field term in Tin is constant 
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for fixed a, the Hamiltonian at a given fixed a is: 



Us = J / dr 



da(r) 



dr 



dr,dr 



cr{ri)c7{rj) 



■ 3 (h-r.l+a)^ ' 
(2) 



limit finite. The function a 1 is periodic with period L, 
and we define its Fourier transform: 



*'(r) - £^ G e 



irG 



°G 



L 



a'(r)e 



-irG 



dr . 



(6) 



(7) 



III. DETAILS OF CALCULATION 

To investigate the T = properties of this classical sys- 
tem, we minimize the energy density to find the ground 
state. We assume that the ground state has a simple 
periodic structure, where each period is comprised of a 
length if of up spins followed by a length if of down spins. 
In the H = case, it has been proved^ that the ground 
state must be of this form, with if = Zj.. In appendix [Al 
Monte Carlo results are presented that confirm a simple 
periodic configuration to be the ground state for H ^ 0. 
A period has total length L = if + if, and simple alge- 
bra reveals that if = (1 + a)L/2 and Zj_ = (1 — a)L/2. 
Minimizing the energy density for a given a under this 
assumption is equivalent to minimizing the energy den- 
sity of a single period with respect to variation in L (we 
choose L, but any of the interdependent variables {if, if, 
L} could be used). The function a(r) is specified by the 
two parameters, a and L: 



cr(r) = 



-1 < (r mod L) < if 
-1 if < (r mod L) < L . 



(3) 



with if as given above. Before writing down an explicit 
formula for the energy density we must choose a zero of 
energy, and the appropriate choice depends on the value 
of A. 



A. Case I. A < 1 

When A < f , we choose the zero of energy to be a uni- 
form spin density of value a. This is equivalent to placing 
the system in a background jellium of "spin charge" , with 
density —a, and results in the replacement of a with 



er'(r) = cr(r) 



(4) 



in |(2J) . The energy density of a configuration with average 
spin a and total period L is given by: 



e(a,L) 



4J 



hm — 

x^oc 4X 



X 



(|r — r'\ + a) 



The first term is the energy of two domain walls per pe- 
riod divided by the period length L. The second term is 
the limit as X goes to infinity, of the energy density due 
to the long-range interaction of a finite system of length 
2X; our zero of energy here has been chosen to make this 



where the sum is over reciprocal lattice vectors G = 
2nm/ L for m S Z. Taking the Fourier transform of the 
second term in ^ gives 



(8) 



G 



We have used that J dr e ir ^ Gl ~ G ^ = 2X8 GuG2 , taken the 
limit X — > oo , and have defined the Fourier transform of 
v(r), vq — f_ dr v(r)e~ lrG . In the case A < 1 we can 
compute vg when v(r) = l/(r + a) x to obtain: 



v G = 2 



cos(Gr) 



o {r + a) 
2G A ~ 1 



dr 



(A<1), (9) 



where we have skipped intermediate steps in the integra- 
tion, and in the second line have let a — > 0. Looking back 
to ijHJ), we next must calculate the Fourier series coeffi- 
cients of cr'(r). We achieve this by calculating the Fourier 
coefficients og of er(r). Then, since a' G = <tg f° r G / 0, 
and <j'g=v> = by definition, we obtain a' G . 



&G 



dr(l)e- 4rG + / dr(-l)e 



-irG 



(10) 



By substituting G — 2-KvnjL, and inserting Q and (flQ]) 
into ([5]). we arrive at a final expression for the energy 
density: 



e(a,L) = j- 



4 J Q2 A+1 r(3- A)sin^ 



L A -!(1-A)(2-A) 
where we have defined 

1 — (—!)"' cos(7rmer) 



C(a,X), (11) 



C(a,X) 



m— 1 



(7rm) 



3-A 



(12) 



which converges for A < 2. Finally, we solve H = to 
obtain the length of the period (L) that minimizes the 
energy density subject to the specified value of a (and 
A < 1) in the limit a — > 0: 



L = (2-A)- 



2 A -!r(3- A)sin^C(a,A) 



(13) 
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B. 



Case II. 1 < A < 2 



When A > 1 it is best to choose the fully polarized 
(ferromagnetic) state as the zero of energy in order to 
eliminate the ultraviolet divergence of the energy inte- 
grals. The energy density for A > 1 then reads: 



x 

e(a,L) — h lim — — f drdr' v (\r — r'\ 

-x 

x (o-(r)o-(r') - 1) . 



(14) 



The first term is unchanged from the previous case, since 
the derivative is only sensitive to domain boundaries, and 
in the second term we explicitly subtract the energy of a 
fully polarized configuration. Performing a Fourier trans- 
form on the second term leads to: 



[a, L) = — + — (v G 



vo) \<tg\ 



(15) 



G#0 



Since oq has been found above (|10p . all that remains is 
to calculate the expression: 



vg 



2G A " 1 r(l - A) sin 



(cos(GV) - 1) 
(r + a) x 
ttA 
2 



dr 



AG [1,2) , (16) 



where the second line, again obtained by taking the limit 
a — > and true only for 1 < A < 2, is precisely the same 
as (0) for A < 1. Since <jg agrees with a' G except at 
G = 0, and G — is excluded in the sum of (fl~5|) . the 
results obtained for A < 1 and 1 < A < 2 separately may 
be combined, expanding the region of validity of (|13p to 
< A < 2. 

In figure[TJ the size of the up spin domain (Z|) is plotted 
at several values of A G (1,2) as a function of a. As 
A increases, the frustrating interaction becomes shorter 
ranged, and the size of the domains becomes larger as 
expected. One must be careful, however, in interpreting 
figure Q] since the unit itself depends on A. This exposes 
a difficulty of working with the continuum model with a 
variable exponent A, as the only length scale is dependent 
both on A and J/Q. Figure [5] shows for fixed A = 3/2 
the size of both spin domains, as well as their sum L, as 
a function of the applied field H (obtained via Legendre 
transform) . 



Case III. A > 2 
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FIG. 1: Size of the spin up domain in the ground state of ([2]) 
as a function of average spin a. The vertical axis is in units 



regions where o(r) ^ o~(r') (i.e. a(r)o~(r') = —1) yield 
nonzero contributions in the second term, and after some 
algebra, for L a the energy density can be written as: 



e(a,L) = 



4J 4Q 
T + (A-2)(A- l)L 



C'{a,\)+a 2 - x 

L \-2 



1 



a A - 2 
(17)" 



where we have defined the following sum, convergent for 
A>1: 



C'(a,\)=J2[(n + a) 



2-X 



2n 



2-A 



+ (n — a) 



2-Al 



(18) 

and a = (1 + c)/2. For larger a/L there are corrections 
to C of order a/L. We set the derivative of the energy 
per site with respect to the period L to zero and thus 
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When A > 2, to keep the energy density finite, the 
ultraviolet cutoff a cannot be set to zero, and a must be 
carried through the calculation. Although the analysis 
leading to (TT5")) and (fTB)) can be extended to A > 2, the 
sum in (|12[) would diverge, and it is preferable to analyze 
the system entirely in real space, starting from (|14p . Only 



FIG. 2: Size of the spin up domain if, spin down domain l±, 
and their sum L (period of the configuration) that minimizes 
the energy of classical Hamiltonian ([2]) with A = 1.5. We 
obtain these quantities as a function of applied field H by 
performing a Legendre transformation. 
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find the ground state L: 

de 4J 4Q 

dZ = "Z 2 " ~ (A-2)L A 
4Q 

(A-1)(A-2)L 2 

-,2-A 

= -J 



,2-A 



A 



J 



1 



A-2 



2 

A- 1 



- L 



2-A 



X - 1 

_ a 2 - x C'{a,X) 
L \-2 



(C"(a,A)+a 2 " A ) 



In the present case of A > 2, C'(a,X) > 0, so the last 
line above will have a solution with finite L only below a 
critical value of J/Q: 



J 



,2-A 



(A-2)(A-1) 



(20) 



For larger values of J/Q the energy is minimized by in- 
finite L so the ground state has macroscopic phase sepa- 
ration, since the energy of a domain wall is then always 
positive. 

Thus, for A > 2, this model with a given J/Q may or 
may not show finite domains in its ground state, depend- 
ing on the cutoff a. In this regime, finite domains form 
in the ground state only for J/Q < (J/Q) c , as specified 
by (|2"0]) . Above this value the ferromagnetic interaction 
is too strong relative to the antiferromagnetic interaction 
for domains to form. Figure [3] shows the regions of phase 
space where domains exist, and where the domain size is 
cutoff independent in the limit a — > 0. 
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FIG. 3: Zero temperature phase diagram of Hamiltonian 1(2). 
showing where finite domains are present in the ground state, 
and whether the domain sizes are dependent on the UV cutoff 
a in the limit a — > 0. 



D. 



Case IV. A 



In the special case A = 2, an analysis similar to the 
above A > 2 case can be done, which results in an energy 



density (valid for L ^ a) 



(- n 4J 



AQ(3 2 



L 

AQa 2 r 



/ a \ aL 
\aL J a 



L 

Qa(3 
L 



, f a \ (3L 



111 



a/3L 



a 
L 



(21) 



(19) where (3 = (1 — er)/2 and F(a) is a function indepen- 
dent of L. By solving = we find that the ground 



state period Lq is exponentially dependent on J/Q with 
a prefactor proportional to a: 



Lq = a exp 



J 

Q 



F'{a) 



(22) 



where F'(a) is another function independent of L and of 
order unity for a G [— 1,+1]. This form can be argued 
from dimensional analysis, since J/Q is dimensionless at 
A = 2, and the only length scale in the problem is a. 
The boundary case A = 2 separates the regime where L 
approaches a finite limit given by (fl~3|) as a — ► for A < 2, 
from the the cutoff-dependent regime for A > 2. 



IV. CONCLUDING REMARKS 

In summary, we have studied the generalized Coulomb- 
frustrated Ising model in a one dimensional continuum 
for different exponents A of the long range interaction 
(A = 1 for the Coulomb case). We have derived an an- 
alytic solution for the ground state domain configura- 
tion of this model under the assumption that the ground 
state configuration has a simple periodic structure. This 
assumption has been confirmed by Monte Carlo simula- 
tions (see Appendix [XJ . Such simulations can be done 
in any dimension d, and mark an avenue for future work; 
more complicated domain patterns presumably do occur 
for d > 1. 

We find that for < A < 2, as the magnetic field H is 
increased from — oo, the ground state is ferromagnetic un- 
til a critical field — H c is reached, at which point isolated 
domains of flipped spins of finite length are formed in an 
otherwise polarized background. For — H c < H < H c , 
periodic configurations with if up spins followed by Z| 
down spins become the ground state, and the system is 
said to be in a "microemulsion" phase. In higher di- 
mensions much numerical work has been done in the ab- 
sence of a magnetic field(see— and applications have 
been made to magnetic thin films 14 in 2D, as well as the 
metal insulator transition in 2D and 3Di^. Analytical 
expressions analogous to our results can be derived as- 
suming the ground state is simply periodic along only one 
direction 8 . But more generally one is presumably forced 
to resort to approximations and numerics. The numeri- 
cal Monte Carlo work, however, can be readily extended 
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to higher dimensions, and can be used to investigate do- 
main formation and behavior as a magnetic field is varied 
(i.e. in the non-charge-neutral case). 

Our solutions describe how Zf increases and Zj de- 
creases with increasing H. At zero field Zf = Zj_, and at 
H = ±H C we find that l-\ or Zj_ diverges to infinity, respec- 
tively, while the length of the "minority" domain remains 
nonzero and finite. For A > 2, whether or not the mi- 
croemulsion phases appear in the transition between fully 
polarized up and down states depends on the dimension- 
less quantity (Qa x ~ 2 )/J. In the case that microemulsion 
phases do occur in this region, their properties depend 
on the ultraviolet cutoff a. Our results provide more 
explicit examples of models with frustrating and suffi- 
ciently long range interactions that have ordered "mi- 
croemulsion" phases instead of macroscopic phase sep- 
aration. However, whether such a result is obtained in 
an inherently quantum system such as the 2D electron 
gas requires further investigation. In this context, it is 
worth noting that claims exist in the literature in favor of 
the classical scenario, both for the electron gas^ and for 
the highly disordered Anderson model with long range 
Coulomb interactions^. 



APPENDIX A: MONTE CARLO ANALYSIS 

The approximation central to this paper is that the 
ground state of TL S (see eq. [2J has a simple periodic 
structure. While this has been proved for H = (i.e. 



a = 0), no such result exists for nonzero magnetization. 
Thus, to justify the approximation, we have performed 
Monte Carlo simulations on a discretized version of eq. 

m 

H 5 = -J^ W+1 + QE3l' ( A1 ) 

i ij 

If we introduce a lattice spacing parameter a, then in 

the limit J/Q -> and a -» such that Qa> ^ 2 = C 
for constant C, the above Hamiltonian is equivalent to 
the continuum formulation (eq. [3J) with J/Q — C. We 
solve for the ground state the discretized model using a 
Monte Carlo algorithm with simulated annealing. Up- 
dates are determined by a Wolff cluster method^ which 
preserves the overall magnetization, and thus a is fixed 
during a simulation run (similar to the analytical cal- 
culation). Enough runs are averaged over so that the 
variation in the sizes of the resulting domains is negli- 
gible. Additionally, the system configuration converged 
upon by each Monte Carlo run is compared to the exactly 
periodic configuration and it is seen that the energy den- 
sity of the exactly periodic configuration is always equal 
to or lower than the energy density of the Monte Carlo 
result. In summary, we find that at any a € [—1, 1], and 
for all values of J / Q the ground state consists of uniform 
spin up domains of length Z} interleaved with uniform 
spin down domains of length Zj_ (but Zj ^ Zj). Detailed 
results will be provided elsewhere^. 
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